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LIST  OF  SYMBOLS 


M-o  mass  of  the  bird,  grams 

Mot     mass  of  the  sleeve,  grams 

M      mass  of  the  spring,  grams 

a      dimensionless  factor  describing  an  effective 
spring  mass 

k      spring  constant,  dynes/cm. 

F      frictional  force  opposing  the  motion  of  the  spring, 
°     dynes 

u      frequency,  radians/sec. 

4a     arithmetic  decrement  of  the  spring  per  cycle,  cm. 

X,      displacement  of  the  bird  mass  with  respect  to  the 
sleeve ,  cm. 

I      represents  an  idealized  impulse  delivered  to  the 
bird  mass  twice  each  cycle,  dynes 


X 

o 


the  distance  of  the  bird  mass  from  equilibrium, 
caused  by  the  rotation  of  the  sleeve  on  the  shaft, 
also  where  the  impulses  occur  (x^  =  ±x^)  ,  cm. 


T      dimensionless  variable  related  to  time  (x  =  ut) 

X      dimensionless  variable  related  to  displacement 
(x  =  x^/x^) 

f      dimensionless  variable  related  to  friction 
°      (f  =  F  /kx  ) 

^  O      0     0^ 

I      dimensionless  variable  related  to  the  impulses 

°      (I   =  I/kx  ) 
^  o     '   o^ 

b      dimensionless  variable  related  to  those  positions 
°     of  the  bird  mass  with  respect  to  the  sleeve  where 
it  locks  and  unlocks  and  where  the  impulses  take 
Dlace . 
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R      radius  of  circular  arc  in  phase  space 

v      dimensionless  speed  in  phase  space 

2 

h      dimensionless  energy  associated  with  energy  given 
to  the  bird  mass  during  an  impulse 

A      the  real  amplitude  of  the  motion  of  the  bird  mass,  cm, 

y      velocity  of  bird  mass  in  phase  space 

Ad     the  amount  the  sleeve  slips  down  the  shaft  each 
half  cycle  as  a  result  of  the  unlocking  of  the 
sleeve ,  cm. 
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CHAPTER  I 
INTRODUCTION 

In  the  celebrated  course  at  M.I.T.  on  nonlinear  mechanics 
J.  P.  Den  Hartog  demonstrates  a  toy  "woodpecker"  to  his  class. 
The  toy  can  be  idealized  as  shown  in  Figure  I.  The  solution 
is  given  as  a  problem  on  the  course's  final  exam  and  Den  Har- 
tog' s  analysis  is  found  in  the  Study  Guide  [1]  for  the  course. 
The  problem  and  his  analysis  are  included  in  all  pertinent  de- 
tail   in   Appendix   A.      His    analysis    concludes   with   the    equation 

f^((}))    =    a^    +    6^    -    Y^    -    5^    =    0  (1-1) 

where 

a  =   A^Ccosh   A(j)-1}  (1-2) 

6    =   A(A(j)-sinh   Ac}))  (1-3) 

Y    =    (A^  +  2)(cosh   A(f)-1)     -    A(})    sinh   Ac})  (1-4) 

6    =    -A   sinh    A<t>   +    A^c})    cosh    Ac})  (1-5) 

with   A   and    c})    defined   in   Appendix   A.      The    solution   of    (1-1) 
would  mean   the    solution   to    the    analysis,   but    this    has   not   been 
done    until    this    point.       In    fact,    Den   Hartog   remarks    [1]    that, 
"This    I   have   not    done,    and   if   anyone    ever   does    it,    I    shall   be 
grateful    if  he    or   she   will    send  me    the    solution." 
Letting 

X   =    cosh   Ac}) 


fv 


(a) 


'SL 


(b) 


Figure    I 


and 

y   =    sinh   A<}) 

and   substituting   these    into    (1-2),    (1-3),    (1-4),    and    (1-5), 
and   thence    into    (1-1)    we   have,    after   canceling    and   gathering 
terms : 

4    7  4    2    2  3  "^ 

f^(()))    =    A^c})^    -    A^4)^x^    -    4A-^4)y   +    AA^<i>xy 

+   A^(-4x^  +  8x-4-(t)^y^)    +   A(A<t>xy-A<i>y)    -    4x^    +    8x    -    4. 

The    right-hand   side    can   be    factored    and   then   both    sides 
divided  by   x-1,    giving: 

^A*^*^  4    2  3  2  cb^v^ 

-^ZT-  =   A%^(-x-l)    +    4A-^(J)y    -A^[4(x-1)    +    |-^] 

+    4A({)y    -    4  (x-1). 
This    expression    can   be    regrouped   and   divided  by    4: 

^A^*^  A%4  ■ 

4^^.1)    =    -x(^  +    A-^-Hl)    +    y(A-^(},+A((.) 

-    4^    (1    -    ^^)    -    a2    .    1. 
(x-l)A^ 

2       2  2 

Noting   that    x    -y      =    1    or   y      =    (x+l)(x-l)    and   dividing   by 

a2    -, 

A  +1  we    get 

,4.2 


^A^*^  A%^       .     ..     .    x.l. 


^ = ^^^!—    [x-Hl    +    ^]    +    A())y    -     (x-1). 

4(x-l)(A^  +  l)  4(A'^  +  1)  A^ 

Rearranging   this    results    in 

^ » =    -    ^^    (x+1)    -H    A<l)y    -    (x-1)    =-    f*(*).  (1-6) 

4(x-l)(A^+l)  ^  ^ 


Clearly  ^ \i 'P)    and  f^((l))  have  the  same  nontrivial  roots 
so  we  need  only  examine  f^(4)),  the  interesting  factor  of 
f^((})).   By  letting  z  =  Ac)) ,  £^(4))  can  be  written  as 

f*((j))  =  z  sinh  z  -  (|)^(cosh  z+1)  -  (cosh  z-1).         (1-7) 

Equation    (1-7)    can   be    explored  by   expanding   sinhz  and   cosh    z 

in    their    familiar  power   series, 

3         ^5  2n+l 

sinh    z    =    z    -H    |j-  ^    3^  -H    .  .  .    -H   ^2lTTTyT  ^    •  •  • 

2  4  2n 

cosh    z    =    l   +    Jy+4T-+...+|^+... 

and  then  substituting  them  into  (1-7),  obtaining 

6      6      6 

^2n+2      ^2n+2     z2(n+l) 
'^(2n+l)  !  "  4(2n!)  "  2(n+l)  !^  "^  

The    first    two    groupings    are    equal    to    zero    and   the    third 
grouping   is   negative    for   all    values    of   z    except    for    z    =    0, 
when    it    is    zero.       The   n        grouping    can  be  written    as 

z2^"2         11  1 


2n!       ^Zn+l    "    4    '    (2n+2 )  (2n+l) ^ 


and  will  be  negative  for  all  n  >_  2  for  all  nonzero  z.   Thus 
f*((J))  is  negative  for  all  z  with  the  exception  of  the  physic- 
ally trivial  case  of  z  =  0,   One  must  conclude  that  the  Den 
Hartog  analysis  is  not  satisfactory  and  the  basic  physical 
model  must  be  reexamined. 


Several  models  of  the  toy  woodpecker  were  built  and 
examined.   The  first  models  built  were  like  that  shown  in 
Figure  I-a.   These  were  not  stable  and  could  not  be  con- 
structed in  such  a  way  as  to  produce  periodic  motion.   While 
this  may  not  be  impossible,  it  seems  quite  difficult,  with 
several  parameters  being  critical.   After  several  hundred 
variations  were  attempted  unsuccessfully,  this  model  was 
abandoned,  and  one  like  that  in  Figure  I-b  was  built.   It  was 
made  using  a  spring  from  a  ball  point  pen,  a  piece  of  broken 
capacitor,  a  scrap  of  wooden  dowel,  and  a  shaft  of  aluminum 
rod,  and  delighted  the  author  by  working  perfectly  the  first 
time  it  was  tried.   By  observing  the  motion  for  various  con- 
figurations several  things  became  apparent.   To  make  sure  of 
what  was  happening,  the  motion  was  videotaped  and  then  the 
tape  was  played  one  "frame"  at  a  time,  allowing  intervals  of 
one-sixtieth  of  a  second  to  be  observed.   The  following 
features  of  the  motion  were  clear: 

1.  The  "bird"  mass  oscillates  up  and  down,  above  and 
below  the  sleeve.   When  this  mass  is  above  or  below  the 
sleeve  by  a  certain  amount,  the  sleeve  is  locked  on  the  shaft 
due  to  friction.   In  between  these  points  the  sleeve  is  un- 
locked and  drops  a  short  distance. 

2.  "Hard"  excitation  is  necessary  to  start  the  motion. 
That  is,  a  certain  threshold  energy  input  was  necessary  to 
start  the  characteristic  slip-stick  motion. 

3.  The  motion  is  periodic  with  no  change  in  frequency 
from  start  to  finish. 


4.  The  sleeve  locks  in  both  possible  orientations. 

5.  The  frequency  of  the  system  is  not  the  same  as  the 
natural  frequency  of  the  bird  mass  on  the  spring  with  the 
sleeve  locked,  but  somewhat  less. 

6.  UTien  the  sleeve  is  locked  and  the  natural  motion 
observed  the  decay  in  the  amplitude  is  not  geometric  but 
appears  to  be  arithmetic. 

7.  The  sleeve  moves  a  short  distance  down  (about  0.15 
centimeter)  twice  each  cycle. 

8.  Running  the  woodpecker  in  a  vacuum  rather  than  in  a 
normal  atmosphere  resulted  in  no  change  in  frequency. 

9.  The  steady  state  motion  is  independent  of  the  initial 
conditions . 

10.  None  of  the  physical  parameters  appear  to  be  critical, 

11.  UTien  the  system  is  started  with  a  very  large  initial 
displacement  the  sleeve  is  observed  to  lock  and  unlock,  but 
not  to  drop  for  several  cycles. 

The  above  observations  indicate  that  this  problem  falls 
in  the  area  of  self-excited  oscillations.   More  specifically, 
it  falls  in  the  area  of  mechanical  clocks.   The  theory 
describing  mechanical  clocks  has  been  worked  out,  most 
notably  by  A.  A.  Andronov  and  coworkers  [2,3]  although  other 
people  [4,5]  have  discussed  it  as  well.   As  Andronov  points 
out  [2,  p.  168],  "A  clock  is  an  oscillating  system  which 
maintains  oscillations  whose  amplitude  is  independent  of  the 
initial  conditions.   To  start  the  clock  a  large  initial  im- 
pulse is  usually  needed.   If  the  initial  impulse  is  too 


small,  then  the  clock  comes  to  rest  again."   He  also  notes 
[2,  p.  168]  that 


Any  clock  mechanism  can  be  broadly  divided  into 
three  parts:   (1)  an  oscillating  system,  for  exam- 
ple, a  pendulum,  a  balancewheel ,  etc.,  (2)  a  source 
of  energy  such  as  a  weight  or  spring,  (3)  a  trigger 
mechanism,  connecting  the  oscillating  structure 
with  the  energy  source.   For  fixed  positions  of  the 
oscillating  system  (referred  to,  for  the  sake  of 
brevity,  as  the  pendulum)  the  trigger  mechanism 
acts  and  energy  is  given  to  the  pendulum  during  a 
short  impulse.   In  a  good  clock  the  impulse  is  of 
very  short  duration.   The  trigger  mechanism  usually 
acts  twice  during  a  period  close  to  the  position 
of  equilibrium  where  the  velocity  is  greatest.   It 
is  important  to  note  that  the  instant  of  time  when 
the  trigger  mechanism  begins  to  act  is  entirely 
determined  by  the  position  of  the  pendulum.   In 
addition,  the  manner  of  its  action  and  the  magni- 
tude of  the  impulse  also  depend  on  the  state  of 
the  pendulum.   Consequently,  all  forces  which 
arise  in  the  mechanism  depend  only  on  the  positions 
and  velocities  of  the  separate  parts  of  the  system, 
and  not  on  the  time.   Thus  the  clock  is  an 
autonomous    system. 

There  can  be  no  question  that  the  toy  woodpecker  is  in 

fact  a  mechanical  clock,  and  can  be  analyzed  with  the  aid  of 

clock  theory.   In  Chapter  II,  the  theory  of  the  woodpecker, 

now  revealed  as  a  cuckoo,  will  be  discussed,  and  in  Chapter 

III  a  comparison  of  the  theory  with  some  experimental  results 

will  be  given.   Chapter  IV  includes  a  summary  of  this  work 

and  its  conclusions. 


CHAPTER    II 
THEORY 

Self-excited   oscillations    form   an    important    class    of 
phenomena   including   some    electronic    oscillators,    buzzers, 
electric  bells,    wind    and   string   musical    instruments,    brake 
"squeal,"    some   kinds    of   aerodynamic    flutter,    and   so    forth. 
These   phenomena   are    characterized   by    definite    oscillations 
which    are   properties    of   a   specific   system   and  not    dependent 
on    the    initial    conditions.       One    of   the    most    interesting 
aspects    of   these    systems    is    that    they    can    generate    a  periodic 
process    from   a  nonperiodic   source. 

Clocks    are    a   distinct    class    of  mechanical    self-excited 
oscillations,    and   consist   of   three   parts:       (1)    an   oscillating 
system,    (2)    an   energy   source,    and    (3)    a   trigger   mechanism   or 
escapement    connecting   the    two.       Clock    theory  has   been  worked 
out    [2-5]    and,    although    the   woodpecker   is    different    in 
several  ways    from  the    clocks    discussed   in   the    literature,    the 
theory    can   be    suitably   extended   to    describe    it. 

The   physical    features    of   the    motion   mentioned   in    Chapter 
I    should  now  be   more    carefully   examined   in    the    context    of 
developing   the    most    descriptive    mathematical   model.      Hard 
excitation    is   necessary   to    start   the    motion.      Soft   excitation 
is    defined   as   when    an    arbitrarily   small    initial    displacement 
is   needed   to   start   the    motion    and  hard  excitation    as   when    an 
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initial  displacement  (or  other  appropriate  energy  input)  must 
be  larger  than  a  certain  amount  to  start  the  motion.   The 
hard  excitation  observed  implies  that  the  friction  present 
is  either  of  the  Coulomb  type  (dry  friction) ,  or  dominated 
by  Coulomb  friction.   A  summary  of  an  argument  by  Andronov 
indicating  that  Coulomb  friction  (compared  to  linear  or 
velocity  dependent  friction)  is  the  more  appropriate  type  for 
clocks  is  given  in  Appendix  B. 

The  actual  motion  of  the  system  is  not  complex  and  can 
be  analyzed  in  a  series  of  several  steps.   Beginning  with  the 
bird  mass  in  the  maximum  upward  position  with  respect  to  the 
sleeve,  and  with  the  sleeve  locked  the  features  of  the  motion 
are  as  follows : 

1.  The  bird  begins  to  descend  and  continues  until  it  is 
a  short  distance  above  the  horizontal. 

2.  At  this  point  the  sleeve  unlocks  and  the  whole  sys- 
tem begins  to  descend. 

3.  The  bird  continues  to  descend  with  respect  to  the 
sleeve  and  when  the  bird  is  a  short  distance  below  the  hori- 
zontal the  sleeve  locks.   At  this  point  from  the  frame  of 
reference  of  the  bird  it  would  appear  that  its  velocity  with 
respect  to  the  sleeve  suddenly  increases.   The  idealization 
of  this  is  an  impulse  applied  to  the  bird  at  the  instant  the 
sleeve  locks. 

4.  The  bird  continues  to  its  downward  extreme  and  then 
starts  back  up. 


un 


an 
an 
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5.  l\Tien   the   bird   is    just   below   the   horizontal   the    sleeve 
locks    and    as    the   bird   continues    up   with    respect   to    the 

sleeve,    the    system   again   begins    to    descend. 

6.  The  bird  continues  up  with  respect  to  the  sleeve 
d  when  it  is  just  above  the  horizontal  the  sleeve  locks 
d   again    this    can  be    idealized   as    an    impulse    applied   to    the 

bird. 

7.  The   bird  continues    to    its    upward  maximum   and   the 
cycle   begins    again. 

The    general    form   of   the    differential   equation    describ- 
ing  the    motion    can  now  be   written    as 

Mx,    +   kx,    =    ±F      +1  (2-1) 

11° 

where  each  of  the  terms  will  be  discussed  below. 

The  mass,  M,  is  equal  to  M„  +  aM   for  that  part  of  the 
motion  when  the  sleeve  is  locked.   The  mass  of  the  bird  is 
Mp.  and  M   is  the  mass  of  the  spring.   The  quantity  a  is  a 
correction  factor  used  to  introduce  an  effective  spring  mass. 
It  would  be  measured  by  locking  the  sleeve  and  then  measuring 
the  frequency  of  the  bird  for  various  values  of  M„.   These 
values  would  then  be  fitted  by  the  curve 

^  =  2T  [m^  +  aMj  ^^'^^ 

using  least  squares  to  get  the  best  value  for  a. 

During  that  part  of  the  motion  when  the  system  is  falling 
M  would  be  given  by 


11 


where  M,-,,  is  the  mass  of  the  sleeve.   For  this  situation  the 
natural  frequency  would  be 

^^    =   M  (2-4) 

with  M  given  by  Equation  (2-3). 

The  spring  constant,  k,  has  the  usual  definition  and 
would  be  measured  by  the  normal  static  means.   The  Hookes 
Law  approximation  is  made  and  must  be  substantiated  by  exper- 
imental results.   Although  the  present  theory  is  based  on 
phase  plane  analysis  with  linear  equations  pieced  together 

on  the  phase  plane,  it  would  be  possible  to  extend  it  to  in- 

3 
elude  a  nonlinear  term  such  as  kx  ,  if  it  were  necessary.   To 

handle  such  a  term  an  approach  such  as  the  "phase-plane  delta" 

method  [6,  p.  244]  would  have  to  be  used. 

The  term  ±F^  is  a  quantity  describing  the  frictional 

losses  within  the  spring.   The  sign  is  chosen  as  opposite  to 

that  of  X,.   It  would  be  measured  by  locking  the  sleeve  and 

measuring  the  arithmetic  decrement  of  the  amplitude  per  cycle, 

4a,    and  noting    [2,   p.    153]    that 

F      =    amw      =    ak .  ( 2  -  5 ) 

o 

As  far  as  the  impulse  term,  I,  is  concerned,  one  of  two 
assumptions  is  usually  made  in  clock  theory:   that  either 
CI)  momentum  or  (2)  energy  is  conserved  during  the  impulse. 
Since  in  the  woodpecker  a  uniform  amount  of  gravitational 


potential   energy   enters    the    system   during   each   cycle,    the 
second   assumption    seems    to  be    the    more    reasonable    one.       It 
is    also    the    one    made    more    frequently   in    clock    theory    [6, 
p.     168]. 

There    is    a   slight    amount    of   rotational   motion   present 
due    to    the    sleeve    locking   in    an    "up"    and   a   "down"   position. 
It    is    possible    to    consider    a   treatment   with    two   variables, 
one    linear    and   one    angular,    to    describe    this    motion.      However, 
because    the    angular   motion    is    relatively   slight,    it    is   pos- 
sible   to    do    an    analysis   with    one    variable   by    approximating 
the    effect    of  the    rotational   motion.       To    do    this,    the    motion 
is    divided   into    three    regions : 

Mx't    +   kx,    =    ±F      +1    for    -Xn    <    x,    <    x„  (2-6) 

llo  (J—       1—       U 

Mx'    +   kx'    =    ±F    ,    x^   =    x^   +    Xq    for   x,    >    Xq  (2-7) 

Mx'    +   kx'    =    ±F    ,    x'    =    Xt    -    x„    for   x,    <    x„  (2-8) 

1  1  °'110  10 

Where    the    impulses    are   pictured   as    occurring    at    x,    =    -^n  '    ^^'^ 
at    this    value    of  x, ,    the   position    of   the   bird   is    changed  by 
X|^  ,    an    amount  which    corrects    x,    for   the    slight    rotational 
motion    of   the    sleeve. 

The  equations  (2-6),  (2-7),  and  (2-8)  can  be  simplified 
by  a  change  of  variable.  First  the  equations  are  divided  by 
kx^ ,    and   then    let 

"^  =  (m)  ' 

and 

T    =    ut  (2-9) 
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so   that 

F  F"  "    ^  (2-10) 

where   now  we   have    a   dimensionless    time,    x.      Also    let 
f      =    'O 


^0 
and 

'0  Xr 


(2-11) 


h^    =    ^  =    1.00,  (2-12) 


^0 
and 

^0  =  k^  '  ^''^'^ 

so  that  now  (2-6),  (2-7),  and  (2-8)  take  the  form 


X  +  X  =  ±£q  +  I    for  -1  i  X  <  +1  (2-14) 

X  +  X  =  ifp  +  bp   for  X  >  1  (2-15) 

X  +  X  =  ±fQ  -  bp   for  X  <  1  (2-16) 

This  motion  can  be  shown  with  a  phase  plane  diagram.   The 
velocity,  y,  of  the  bird  with  respect  to  the  sleeve  is  repre- 
sented by  the  vertical  axis  with  upwards  (positive)  y  corre- 
sponding to  downward  motion  of  the  bird.   The  x  axis  repre- 
sents the  displacement  of  the  bird  with  respect  to  the  sleeve 
which  is  also  positive  for  downward  motion  of  the  bird. 
There  are  two  impulses  per  cycle  which  in  this  model  are  con- 
sidered to  be  of  equal  magnitude,  occur  symmetrically  in  each 
half-cycle  and  take  no  time.   The  phase  diagram  is  shown  in 
Figure  II.   The  convention  used  in  clock  theory  is  to  follow 
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Figure    II 
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the  phase  path  in  a  clockwise  sense.   At  point  S  the  bird  is 
at  its  maximum  negative  (upwards)  displacement  and  its 
velocity  is  zero.   The  bird  then  begins  to  speed  up  as  it 
falls.   At  X  equal  to  minus  one  the  sleeve  unlocks  and  the 
system  falls,  the  bird  going  slower  with  respect  to  the 
sleeve  than  before.   At  x  equal  to  plus  one  the  sleeve  locks, 
an  impulse  is  delivered  to  the  bird  and  it  continues  to  its 
positive  (downward)  extremum  at  f .   This  process  repeats  it- 
self in  the  other  direction  ending  at  S.   The  phase  paths  in 
the  various  parts  ot  the  motion  are  portions  of  circles 
centered  at  (0  ,±fr.)    or  (0,±fQ+bp)).   This  can  be  seen  by 
letting  y  =  x;  then  equation  (2-14),  for  example,  can  be 
written  as 

y  +  X  =  fp  (x  7^  ±1) 

or 

V  ^  ^  ^  =  fo- 

Rearranging,   we   have 

ydy   =    (fQ-x)dx 

1         2      2 
and   letting   j   (R   -f^)    be    a   constant    of   integration,   we    get 

after   integration 

y2    =    2fQX    -    x^    +    r2    -    £2 

or 

y2    +    (x-fQ)^=   R^  (2-17) 

which  is  a  circle  in  the  phase  plane  centered  at  (Oj-f^). 
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For  periodic  motion  we  require  that  the  speed  at  point  A 
be  the  same  as  that  at  point  B: 

^A  =  ^B  (2-18) 

This    condition  will    allow   the    frequency    and    amplitude    of   the 
motion    to  be    calculated.      As    shown    in    Figure    III,   we   have 

R-    +    2f      =    R,  ,  (2-20) 

2  o  1  ' 

and 

R^    -    (1-f   )2    =    R^    -    (b    +f    -1)2    .  (2-21) 

3o2oo 

Also,   we    see    that 

(v;)2    +    (f^+l)^    =    4  (2-22) 

and 

(vl)2    +   h^    =    v^  (2-23) 


where 

A(y2)    =   h^  (2-24) 

is    the    change    in   the    square    of   the    speed   in    dimensionless 
units    associated  with   each    impulse.      The    equations    above   must 
be    solved  to    obtain    a   relation  between   v,    and  v„.      Substitutinj 
first    (2-22)    and   then    (2-21)    into    (2-23)    we    obtain 

v^    =  h^    -    (1+f   )2    +    r2    +    (i-f  )    -    (b    +f   -1)2 
B  ^o'^  2^o'^^oo-^ 

and  then  substituting  (2-20)  into  this,  we  have 
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—  X 


Figure    III 
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B  ^o  ^o^  ^00^  1  ol  o 

Using    (2-19)    for   R^    one    gets 

-    Cb„*f„-1)2    *    4(b^*f„)2    -    4(b^*f„)Rj. 

Setting 

2  2  2 

^B   =    ^A   =    ^    ' 

and  canceling,   we    get 

4f  R,    =    4f^    +   h^    -    4b    £    .  (2-25) 

0     1  O  0     0  ^  ^ 

Dividing  by    4£    ,    squaring  both   sides,    and   again   using    (2-20) 

results    in 

2 

v^    +    (1+b    +£    )2    =    £2    -    2£  b      +   ^  +    b^    -    =^  +    J^   . 
^0      0-^  o  oo2  o2£  i/:r2 

o  lot 

o 

Rearranging   terms    gives 

,  ,4  ,2        b   h^ 

v^    =   -^  +   ?-  -  ^ —  +    2b      -    1    -    2£      .                                        (2-26) 

,^r2         2  2£                  o                      0 

16 1  o 
o 

It  is  still  necessary  to  determine  h  to  be  able  to  cal- 
culate V.   By  noting  that  the  system  drops  a  specific  dis- 
tance. Ad,  once  each  half  cycle, the  energy  input  during  the 
impulse  can  be  calculated  by  assuming  it  is  due  to  this 
change  in  gravitational  potential  energy: 

g(Ad)  =  i  (Av)2  (2-27) 


19 


where 


CM^  =  (X?) 


rz 


r-2, 
(x^)^ 


is  the  change  in  the  square  o£  the  physical  speed  as  a  result 

of  the  impulse.   This  must  still  be  related  to  the  dimension- 

2 
less  quantity  h  .   By  definition 


X. 


3T  '-x  h 


X 


^  ^^h 


and  using   equation    (2-9)    we   have 


1    d       .^l. 

oi  W   ^x    ■'2 
o 


1  ^    (% 

CO    It     ^x    -^ 


Then   substituting   in   the    definition   of   u,   we    get 


h"-   = 


M 


kx 


2    [^^1^2 


Cxp,] 


or 


kx 


M  r  /- A      ^  2t 

2"   [(^v)    ] 


and   using    (2-26),   we    finally    get 


h^    =    -^   [2g(Ad)] 


(2-28) 


kx 


o 


The  mass,  M,  in  this  equation  is  M„  +  aM  and  does  not  in- 
clude the  sleeve  mass  since  the  energy  the  sleeve  picks  up 
will  be  lost  when  it  is  stopped  by  the  friction  between  it 
and  the  shaft. 
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Now  that  h  has  been  obtained  it  is  possible  to  solve  for 
V  and  then  to  solve  for  R^ ,  R2 ,  and  R^  using  the  equations 
developed  above.   The  frequency  of  the  motion  can  be  calcu- 
lated by  noting  "that  the  representative  points  move  along 
the  phase  paths  with  an  angular  velocity  equal  to  unity" 
[2,  p.  197].   Thus  the  time  of  transit  in  x  units  along  one 
of  the  semicircular  arcs  is  equal  to  the  angle  subtended  at 
its  center  by  the  arc.   As  shown  in  Figure  III,  the  relevant 
angles  are  a,  3,  6,  and  y,  and  the  period  in  t  units  would 
be  given  by 

T  =  2(a  +  6  +  5  +  Y  +  j)  .  ■  (2-29) 

By  inspection  of  Figure  III  the  angles  above  can  be  easily 
obtained  in  terms  of  known  quantities: 

tan   a   =    ^.^    .f  ^    ■    ^ 

o      o 

b    +f    -1  ,  , 

o      o  f2-31) 

sm   B    =   D ^        ^ 

^2 

sin   6    =   ^  (2-32) 

^3 

1  +  f 

o  r2-33) 

sm   Y    =    -R —  ^         ^ 

*^3 
In   physical    units    the   period  would  be 

and  the    frequency  would  be    1/T. 
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The  energy  balance  is  an  important  aspect  of  this  motion 
For  regular,  periodic  motion  the  energy  input  per  cycle  must 
equal  the  losses  per  cycle.   The  energy  input  for  the  bird 
and  spring  is  on  the  order  of  2 (M^+  M  )g(Ad)  per  cycle  and 
the  energy  losses  per  cycle  are  on  the  order  of  4f  A  or  4akA 
where  A  is  the  amplitude  of  the  motion  measured  from  the 
equilibrium  point.   Equating  these  two  energies  and  solving 
for  the  amplitude,  we  have 

A  =  --fi^  (2-34) 

To  the  extent  that  the  energy  balance  is  correct,  this  rela- 
tion should  predict  the  amplitude.   This  balance  is  basic  to 
the  theory  and  (2-33)  can  be  derived  another  way  by  noting 
that  the  normalized  amplitude,  A  ,  where 

A  =  A/x  , 
o      o ' 

can  be  derived  from  the  phase  diagram,  Figure  III,  as 

A  =  R^  +  b   -  f  . 
o    1    o    o 

Using    (2-25)    for   R,    we   have 

A   =    ^  =    tL 
^o        ^^0 

Substituting  (2-5),  (2-11),  and  (2-28)  into  this  results  in 


A_  ^  Mg(Ad) 

X     2akx 
o        o 


which  is  the  same  as  (2-34) 


Before  turning  to  an  experimental  exploration  of  this 
theory,  some  general  points  about  it  can  be  made.   The  theory 
predicts  that  the  frequency  will  be  less  than  the  natural 
frequency  of  the  bird  on  the  spring  with  the  sleeve  locked. 
We  can  understand  why  this  is  so  physically  because  the  peak 
to  peak  amplitude  of  the  woodpecker  will  be  4b   longer  than 
the  amplitude  of  the  natural  system,  other  things  being  equal, 
and  this  will  take  longer. 

We  also  note  the  necessity  for  at  least  one  dynamical 
measurement.  Ad  being  the  most  appropriate.   In  this  it  is 
like  the  Den  Hartog  theory  where  f,  the  coefficient  of  fric- 
tion between  the  sleeve  and  the  shaft,  is  idealized  as  having 
the  same  value  for  motion  as  for  rest.   In  fact,  this  would 
be  unlikely  and  some  sort  of  dynamic  measurement  would  be 
necessary  to  establish  some  appropriate  average  value  for  f. 

In  the  present  theory  it  is  possible  to  get  around  this 
necessity  for  a  dynamic  measurement,  at  least  theoretically, 
although  some  approximations  will  be  necessary.   First 
assume  that  when  the  sleeve  unlocks  it  falls  freely  with  an 
acceleration  of  g.   Since  it  starts  from  rest  we  have 

Ad  =  i  gt^  (2-35) 

where  t  will  correspond  to  that  part  of  the  period  when  the 
system  is  unlocked: 

1/2 


t  =  C6+Y) 


(    M  Mg^ 
k(M+Mg^) 


(2-36) 
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Now 


6    =    sm         --n 


and 


y   =    sm 


-1  0 


R. 


and   these    can  be    expanded   in    a  power   series 


3  5 

.     -1  X      ^    X 

sm        X   =   X    -    -YY  "*"   TT 


By  taking  only  the  first  two  terms  of  the  series  (this  and 
other  approximations  remain  to  be  justified  in  Chapter  III) 
we  have 


1-f    1+f 


6  +  Y  = 


and 


R. 


'  -^  y  =  k 


R. 


2+6r 


6r: 


(i-f  )^   (i+f  )^ 

^    o         o 


6r: 


6r: 


It  is  further  assumed  that  f   is  small  and  terms  of  the  order 

o 

2 
f   and  higher  can  be  neglected.   So  we  have 


6  +  Y  = 


R, 


3R 


I  • 


(2-37) 


By   substituting    (2-20)    and    (2-19)    and    then    (2-26)    into    (2-21) 
we    obtain, after   some    algebra 


h   ~   W 


o 


8f 


o 


16f^        16f^         32f^]l/2 
o        o        o 


(2-38) 
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Again,    terms    of  the    order    f      and   above    are   neglected,    result- 
ing   in 


8f 


^3   "   3T 


(1 


0)1/2 


o  h" 

Expanding  this  in  a  series  and  keeping  only  the  first  two 
terms  gives  approximately 


(2-39) 


^2       4f^ 

^3  "  ^T~  f^^  "  7T~'^ 
o      h 


or 


R 


3   Tf 


1. 


(2-40) 


Substituting  this  into  (2-37)  yields 
6  +  Y  = 


o     16 


~2  3     2      J 

h^-4f     ^      (h  -4f  ) 
o      ^     c^ 


The  second  term  is  neglected  and  substituting  this  for 
6  +  Y   in  (2-36)  gives 


Ad=f 


M  M 


SL 


k(M+Mgj^) 


1/2  f  8f    "2 
'       o 


h^-4f 


(2-41) 


Letting 

''  M  M 


^    2 


and 


SL 


1/2 


k(M+M^,) 


h^  =  B(Ad) 


where 
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and  putting  these  into  (2-41)  gives 

7    6  4L£^ 

(BAd-4f^)^  =  -^   .  (2-42) 

This  relation  can  be  solved  numerically  to  give  a  Ad  for  a 
particular  M„.   In  Chapter  III  this  equation  will  be  explored 
with  specific  data  to  determine  its  validity. 

Although  only  the  steady  state  motion  has  been  explored 
here  it  would  not  be  difficult  to  extend  the  phase  plane 
analysis  to  cover  various  types  of  transients  present  when 
the  motion  begins.   Several  references  (see,  for  example 
[6,  p.  274,196],  [2,  p. 178],  and  [3,  p.  13])  discuss  various 
types  of  transient  analysis. 


CHAPTER  III 
EXPERIMENTAL  RESULTS 

The  present  work  was  originally  intended  to  be  entirely 
theoretical  in  nature,  consisting  of  completing  the  analysis 
begun  by  Den  Hartog  and  then  exploring  various  physical  situ- 
ations with  the  completed  theory.   Because  it  was  necessary  to 
redo  the  theory,  it  was  felt  to  be  prudent  to  include  experi- 
mental results  with  the  new  theory  so  that  a  complete  and  con- 
sistent picture  would  be  present  for  the  first  time. 

Because  the  original  thrust  of  this  work  was  not  experi- 
mental, the  measurements  were  made  in  most  cases  on  the  basis 
of  rather  limited  resources.   Fortunately,  most  of  the  physi- 
cal quantities  required  are  not  difficult  to  measure  and  the 
results  form  a  consistent  pattern  which  unquestionably  sub- 
stantiates the  theory. 

The  frequencies  of  the  woodpecker  model  for  various  bird 
masses  were  measured  first.   They  were  measured  by  counting 
forty  or  fifty  cycles  and  timing  them  with  a  stopwatch.   Many 
runs  were  made  with  the  same  bird  mass  with  different  numbers 
of  cycles  (for  example,  10,  20,  40  and  50  cycles)  to  see  if 
the  woodpecker  slowed  down  or  speeded  up  as  it  moved  down  the 
shaft.   No  such  effect  was  observed.   These  measurements  as 
well  as  many  others  were  also  made  through  the  use  of  a  video 
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tape  arrangement  made  available  for  this  work  from  the  Orlando 
Division  of  the  Martin-Marietta  Corporation.   The  woodpecker 
was  taped  in  front  of  a  numbered  centimenter  grid  with  a 
digital  clock  reading  in  hundredths  of  a  second  also  in  the 
picture,   because  the  video  tape  could  be  played  back  and 
"framed"  to  show  a  still  picture  of  the  motion  every  sixtieth 
of  a  second.   Appendix  C  contains  some  examples  of  data  col- 
lected from  the  video  tape. 

Table  I  shows  the  average  frequency,  f  ,  obtained  for 
various  bird  masses,  M„.   The  masses  were  measured  to  an 
accuracy  of  about  0,05  gram.   The  individual  frequency  measure- 
ments were  only  precise  to  about  0.1  Hertz,  but  because  large 
numbers  of  measurements  were  made  the  error  estimate  was  cut 
in  half. 

Table    I 

Mg   ±    0.05    (grams)  f^    ±    0.05    (Hertz) 

0.00  4.32 

10.15  2.95 

11.50  2.83 

12.00  2.77 

13.50  2.70 

16.20  2.59 

18.10  2.49 

20.70  2.42 

21.20  2.32 

24.55  2.21 


The  spring  constant,  k,  was  measured  in  the  usual  way, 
measuring  deflection  for  various  loads.   There  was  some  scat 
ter  and  a  slight  indication  of  spring  "hardening"  at  large 
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deflections.      Table    II    gives    the   values    of   the    loads    in   grams 
and   the    respective    deflections. 


Table    II 

Load  ±    0.05 

Deflection   ±    0.05 

(grams) 

(centimeters) 

2.00 

0.40 

3.00 

0.60 

4.00 

0.75 

5.00 

0.85 

5.93 

1.10 

6.00 

1.10 

7.00 

1.25 

8.00 

1.45 

8.43 

1.50 

9.00 

1.50 

9.68 

1.68 

10.00 

1.70 

10.15 

1.90 

10.90 

1.87 

11.00 

1.75 

11.37 

2.10 

12.19 

2.04 

13.00 

2.10 

15.00 

2.30 

17.71 

3.01 

18.96 

3.22 

20.18 

3.44 

21.45 

3.54 

The    average   k    for   these   values    is    5.738   grams    per    centimeter 
or   5,623    dynes    per   centimeter.      A  graph    of   the    loading   versus 
deflection    is    shown    in   Figure    IV. 

The    effective    spring  mass    coefficient,    a,    could  now  be 
determined  using   the   value    for   k   obtained   above.      This    quan- 
tity was    determined  by   locking   the    sleeve    in  place    and  measur- 
ing  the    frequency,    f,    for  various   M„,      Then    a  was    determined 
from  the    equation 
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t: 


-7 


g 


DEfLECTIOKl     IN   Ch. 


Figure    IV 
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f  =  1 


2TT  ^M  +aM 


1/2 


or,  solving  for  a: 

47T  f  M      S 

The  mass  of  the  spring,  M„ ,  is  5.1  grains  and  the  following 
values  for  a  were  obtained: 


Mg  ±  0.05 

f  ±  0.01 

a 

(grams) 

(Hertz) 
3.19 

9.70 

0.84 

11.30 

3.03 

0.83 

12.  75 

2.86 

0.91 

17.75 

2.54 

0.85 

18.99 

2.46 

0.90 

The  average  value  of  a  is  0.87  and  this  was  the  value  used. 

The  quantity  Ad  was  measured  for  each  of  a  variety  of 
bird  masses  by  looking  at  the  total  distance  the  system 
dropped  for  a  certain  number  of  cycles  and  then  averaging  to 
find  the  distance  per  half  cycle.   By  inspection  of  the  video 
tape  it  was  determined  that  the  drops  were  uniform  although 
the  drops  when  the  bird  mass  was  going  down  were  slightly 
different  from  those  when  it  was  going  up.   Table  III  summar- 
izes the  values  of  Ad  for  the  respective  bird  masses. 

The  most  difficult  quantity  to  measure  accurately  was  4a, 
the  arithmetic  decrement.   This  quantity  was  measured  by  lock- 
ing the  sleeve  in  place,  starting  the  bird  mass  oscillating, 
and  then  measuring  the  amplitude  from  the  equilibrium  point 
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Table 

III 

Mb    ±    0.0  5 

Ad 

±    0.005 

(grams) 

(centimeters) 

0.00 

0.2  31 

10.15 

0.182 

11.50 

0.175 

12.00 

0.173 

13.50 

0.166 

16.20 

0.152 

18.10 

0.143 

20.70 

0.130 

21.20 

0.128 

24.55 

0.112 

for   each   half   cycle.       During   these   measurements    the   bird  mass 
would  begin  precessing  which   made    the    vertical    amplitude 
appear   to   decay  more   quickly   on   the   video   tape    than   it    actu- 
ally was.      To    avoid   as    much    error   as   possible,    only   the    first 
five    or   six   extrema  were    considered.      For   a  bird  mass    of    12.75 
grams,    for   example,    the    following  values  were    observed: 

Run    I  ■ 

Extrema     Amplitude  2a 

(centimeters)    (centimeters) 


up 

5.179 

down 

5.083 

-0.096 

up 

4.986 

-0.097 

down 

4.839 

-0.147 

up 

4.451 

-0.388 

down 

4.562 

+0.111 

up 

4.104 

-0.458 

down 

3.9  39 

Run    II 

-0.165 

down 

6.12 

up 

5.93 

-0.19 

down 

5.98 

+  0.05 

up 

5.83 

-0.15 

down 

5.75 

-0.08 

up 

5.32 

-0.43 

down 

5.59 

+  0.29 

Of  course,  many  more  runs  were  made  and  analyzed  but  these 
data  are  shown  as  examples  of  the  data  obtained,  showing  the 
experimental  scatter  and  the  precessional  effects.   It  was 
not  possible  to  make  a  defensible  statement  with  regard  to 
the  nossible  variation  of  a  with  M„,   Since  it  was  not 
desired  to  fit  a  to  theoretical  data,  a  value  of  0.044  centi- 
meters was  assigned  to  all  values  of  M„. 

It  is  possible  that  the  damping  could  be  linear  (propor- 
tional to  velocity  in  some  way)  rather  than  constant.   The 
woodpecker  system  was  run  in  a  partial  vacuum  (<0.1  atm)  to 
see  if  there  was  any  velocity  dependent  damping  due  to  air 
resistance.   No  effect  was  observed  and  the  system  ran  exactly 
the  same  (frequency  and  amplitude)  as  it  did  in  a  normal 
environment.   It  is  still  possible  that  the  losses  in  the 
spring  could  be  due  to  linear  friction  and  this  should  be 
explored. 

In  the  case  of  linear  friction  the  ratio  of  two  succes- 
sive maxima  would  be  given  [6,  p.  201]  by 

n+1     -27TV  ,_  _, 

— —  =  e  (3-2) 

n 
or 

^  =  I¥  Ln  5^  (3-3) 

n+1 

where  v  is  called  the  damping  ratio,  and  where  the  damping  is 

small.   In  the  case  of  constant  friction  we  have 


X   -  ^^,  +  1/2    =  2a  =  constant.  (3-4) 
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Substituting   some    experimental   values    into    (3-3)    we   have 

1      ,       5.18 
"^   =    2lr   L"    4799" 

or   V   =    0.0059    to    the    correct    order   of  magnitude.      With   this 

—  O  7T\)  O 

value    of   V,    e  is    0.964.       If  we    expand   e'    ^^   in    a  power 

series,    consider   only   the    first   two    terms,    and   substitute 
in    the    argument  we   have 

e"2^^   =    1-0.037   =    0.963  . 


Hence,  to  a  good  approximation 

%^=  1  -  V  (3-5) 

n 

and 

X   -  X  ^,  =  vx  (3-6) 

n    n+1     n  ^        ■' 

for  a  full  cycle  and 

^n  "  ^n  +  1/2  ^    2    \  ^^''^^ 

for  a  half  cycle.   In  the  woodpecker  we  note  that  at  the  end 
of  each  half  cycle  the  impulse  restores  the  amplitude.   In 
fact,  it  is  this  exact  relationship  of  energy  balance  which 
gives  the  woodpecker  its  clocklike  behavior.   We  note  the 
resemblance  of  (3-7)  and  (3-4)  for  small  v  and  the  wood- 
pecker's physical  situation  and  conclude  that  whatever  linear 
friction  is  present,  if  any,  can  be  treated  as  constant  to  a 
good  approximation. 

The  final  experimental  quantity  necessary  to  be  measured 
is  x  ,  the  displacement  of  the  bird  mass  with  respect  to  the 
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horizontal  at  which  the  sleeve  locks  and  unlocks  and  at  which 

the  impulses  take  place.   This  was  measured  statically  by 

rotating  the  sleeve  up  and  down  to  its  two  locked  positions 

without  letting  the  sleeve  move  vertically  and  measuring  the 

resulting  displacement  of  the  bird  mass.   These  measurements 

give  an  x   of  0.75  ±  0.01  centimeter. 
*        o 

The  final  results  of  the  frequency  calculations,  a  sum- 
mary of  data  values  used  to  get  them,  and  the  experimental 
frequencies  are  shown  in  Table  IV. 

Table  IV 

k  =  5,623  dynes/cm    Xq  =  0.75  cm 

a  =  0.044  cm  ctMsL  =  4.43  grams 


Mg    (grams) 

d    (cm) 

f 

f 
— e- 

%    error 

0.00 

0.231 

4.50 

4.32 

-4.1 

10.15 

0.182 

2.90 

2.95 

+  1.6 

11.50 

0.175 

2.85 

2.83 

-0.7 

12.00 

0.173 

2.81 

2.77 

-1.4 

13.50 

0.166 

2.70 

2.70 

0.0 

16.20 

0.152 

2.53 

2.59 

+  2.4 

18.10 

0.143 

2.43 

2.49 

+  2.4 

20.70 

0.130 

2.  30 

2.42 

+  4.9 

21.20 

0.128 

2.28 

2.32 

+  1.7 

24.55 

0.112 

2.15 

2.21 

+  2.8 

The  agreement  seems  to  be  very  good.   By  considering 
the  values  of  the  amplitude  predicted  from  (2-34), 

_  Mg(Ad) 

we  can  understand  better  one  source  of  error  in  the  frequency 
calculations.   The  quantity  a  is  the  most  inaccurate  in  (2-34) 
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and  is  most  responsible  for  error.   Table  V  shows  the  calcu- 
lated amplitudes  compared  with  experimentally  observed  ones. 
The  experimental  amplitudes  here  are  taken  from  the  equilibrium 
position  and  obtained  from  a  least  squares  fit  to  five  values. 
The  agreement  is  not  bad,  although  it  could  be  improved,  most 
probably  by  using  more  accurate  values  of  a. 


Table 

V 

Mg    (grams) 

Theoretical 

Exp 

erimental 

%    error 

A    (cm) 

A    (cm) 

10.15 

5.26 

5.40 

2.6 

11.50 

5.52 

5.49 

0.5 

12.00 

5.63 

5.52 

2.0 

13.50 

5.89 

5.62 

4.9 

16.20 

6.21 

5.80 

7.1 

18.10 

6.38 

5.93 

7.6 

20.  70 

6.47 

6.11 

5.9 

21.20 

6.50 

6.14 

5.8 

24.55 

6.43 

6.36 

1.1 

The  last  step  in  the  calculations  will  be  to  examine  the 
validity  of  the  approximate  formula  for  Ad  derived  at  the  end 
of  Chapter  II : 

.2 


64cf' 


(BAd-4f^) 


o 


"AT 


(2-42) 


where 


and 


c  =  g 
^    2 


M  M 


SL 


k(Mg^+M) 


1/2 
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Calculating   Ad   for    an   M      of   20.70    grams    one    gets    a  value 
of   0.23   centimeters,   which    is   higher   than   the    observed   value, 
and  with    this    value    of   0.23   a    frequency    of   2.34   Hertz.      This 
is   higher   than   the    original    calculated   frequency  but    is    in 
better    agreement  with    the    observed    frequency   of   2.42    Hertz. 

Each    step    of   the    derivation    of    (2-42)    can  be    examined   to 
see   how    appropriate    the    approximations    are.      With   the    original 
data   for   an  M„    of   20.70    grams,    the    actual   value    of   R,    is 
found   to  be    7.626.      The    approximate    formula    (2-40), 

h^ 
o 

gives  a  value  of  7.622  which  is  in  excellent  agreement.   The 
actual  value  of  6+y  is  found  to  be  0.263,  whereas  the  approxi- 
mate formula 

6  +  Y  = 


h^-4f 
o 


gives    a  value    of   0.262   which    is    also   in  very   good   agreement. 
The    main   source    of  error   responsible    for   the    two   values    of  Ad 
not    agreeing    any  better   is   probably    in 

Ad   =    i  gt^ 

in  that  the  system  is  not  in  free  fall  but  accelerating  at 
less  than  g  due  to  the  presence  of  some  friction  between  the 
sleeve  and  the  shaft  even  when  the  sleeve  is  unlocked.   A 
smaller  value  of  the  acceleration  would  make  the  calculated 
value  of  Ad  smaller  which  is  the  correct  direction.   An 
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accurate  measurement  of  the  actual  friction  between  the  shaft 
and  the  sleeve  would  call  for  rather  sophisticated  dynamic 
measurements  and  is  probably  not  worth  doing  in  the  sense  that 
Ad  can  be  measured  directly  much  more  easily.   In  fact,  it  is 
a  characteristic  of  this  type  of  problem  that  the  friction 
present  is  idiosyncratic  and  complex.   It  must  usually  be 
idealized  to  a  large  extent  for  a  class  of  problems  rather 
than  dealt  with  for  a  specific  case.   In  this  sense  the  analy- 
sis here  seems  to  be  one  of  the  few  cases  where  a  real  situa- 
tion involving  impact  maintained  oscillations  is  studied 
completely. 


CHAPTER  IV 
COMMENTS 


This  work  begins  with  an  analysis  by  J.  P.  Den  Hartog 
and  shows  that  this  approach  to  solving  the  problem  o£  the 
"toy  woodpecker"  is  not  satisfactory.   Several  possible  reasons 
for  the  failure  of  this  analysis  are  discussed  although  a  com- 
plete "negative"  analysis  indicating  the  exact  reasons  for 
the  failure  of  the  model  and  proposed  solution  has  not  been 
done.   This  work  proposes  a  different  physical  model  for  the 
woodpecker  than  that  chosen  by  Den  Hartog  and  one  which  corre- 
sponds more  closely  to  the  toys  actually  observed.   This 
physical  model  is  then  analyzed  mathematically  and  the  fre- 
quency and  amplitude  of  the  motion  predicted  for  various 
physical  situations. 

A  major  step  in  the  theoretical  development  is  to  recog- 
nize the  toy  woodpecker  as  belonging  to  a  class  of  phenomena 
known  as  mechanical  clocks.   Clock  theory  has  been  worked  out, 
most  notably  by  the  Russian,  A.  A.  Andronov,  and  his  coworkers. 
This  theory  was  extended  and  applied  to  the  toy  woodpecker. 
The  method  involves  piecewise  linear  differential  equations 
matched  together  on  the  phase  plane  and  is  very  suitable  for 
handling  the  impacts  assumed  in  the  model. 
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Actual  working  models  of  the  toy  were  built  and  measure- 
ments of  several  variables  of  the  models  were  taken  and  com- 
pared with  the  theoretical  results.   The  agreement  is  very 
good  and  it  is  clear  that  the  theory  successfully  describes 
the  physical  motion  of  the  toy  woodpecker.   The  calculations 
and  measurements  are  not  difficult  to  do  for  the  most  part 
and  it  is  suggested  that  the  toy  woodpecker  model,  appropri- 
ately constructed,  would  make  a  useful  mechanical  analogue 
for  dealing  empirically  with  various  types  of  self -excited, 
impact  maintained  oscillations  described  by  equations  like 

(2-1). 

Although  the  analysis  is  necessarily  more  complex,  the 
theory  can  be  extended  to  cover  various  other  situations  with- 
out too  much  difficulty.   If,  for  example,  we  consider  a  situ- 
ation with  a  nonlinear  spring  we  could  have 

Mx  +  k^x  +  k2X^  =  ±F  .  (4-1). 

This  could  be  analyzed  in  the  same  fashion  as  that  in  Chapter 
II  if  some  changes  are  made.   The  arithmetic  decrement  must 
be  recomputed  for  example.   This  can  easily  be  done.   The 
potential  energy  for  a  spring  with  the  nonlinear  spring  con- 
stant is 

,2,4 

V  =  -i—  ^  -T-  ^"^'^^ 

The  work,  W,  spent  in  overcoming  friction  will  not  depend  on 
the  velocity  (for  Coulomb  friction)  but  only  on  the  path. 
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If  X    is  the  amplitude,  measured  from  equilibrium,  of  the 
ol 

first  maxima  and  x  ^    the  amplitude  of  the  first  minima,  the 
work  done  in  overcoming  friction  will  be 

The  potential  energies  at  x^^  and  Xq2  are 

2    ,4 
,-   _  ^1^01  ^  ^2^01 
^1 2        4 

and 

,2    ,4 
^1^02  ^  ^2^02 


2      2        2 
from  (4-2).   By  conservation  of  energy  we  know  that 

or,  substituting  in  we  have 

— 2 — 4 2 4      U^oi'  1^0  2''"^  • 

For  k   equal  to  zero  we  have  the  situation  for  the  linear  case, 
I    II    I    2F    ^^  r4-41 

with 

2a   =  Ix   I  -Ix   I  (^4-5) 

^^L  ~  '  Ol'  '  02' 

as  noted  in  Chapter  II.   For  \^    not  equal  to  zero  we  have 


^  (x^l-^02^  '   "t-   ^AvAl^    -    U^oil^l^Oll^P 
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Canceling  like  factors  results  in: 
k  k 

This  can  be  factored  and  written  as 

^2   ,1. .   ,2.  ,    ,  2^-,    2F 


fl^oil-l-ozl^ti  *  2I7  Uxoil'^lxo^r)]  =  i^ 


Substituting  in  (4-5)  and  defining 


^  F 
^N.L.  "  k^ 


(4-7) 


as  the  nonlinear  value  for  the  arithmetic  decrement,  we  get 

^N.L.    -    ^[1  *    W^   (Uoil'*|x„2|2)]  (4-8) 

as  the  formula  relating  a  to  the  observed  amplitudes. 
It  is  possible  to  solve  (4-1)  on  the  phase  plane 
although  a  more  involved  process  of  geometrical  construction 
is  required.   First,  the  equation  is  rewritten  in  the  form 

1       7   ■^  -  F 
and  letting 

we  have 

••^r2   ^^2   3-F     2.^  ^^   -    n 


or 

X  +  p^(x+6)  =  0  (4-9) 


where 

P 

Equation  (4-9)  is  now  in  the  standard  form  for  solution  by 
the  phase-plane  delta  method  [6,  p.  244].   This  method  in- 
volves approximating  6  as  a  constant  over  short  intervals  of 
the  motion  and  constructing  the  phase  paths  as  short  arcs 
approximating  the  true  path. 

Linear  friction  can  also  be  included  in  a  very  straight- 
forward way  if  it  is  proportional  to  the  simple  velocity 
[2,  p.  170].   For  linear  damping  of  a  type  proportional  to 
X  where  n  is  not  equal  to  one  and  may  be  fractional,  more 
sophisticated  methods  may  be  useful  although  the  phase-plane 
delta  method  can  still  be  used.   For  small  nonlinearities  a 
variety  of  methods  are  available  [3,4]  and  can  be  applied  as 
necessary.   On  the  whole,  it  is  remarkable  that  the  toy  wood- 
pecker can  be  dealt  with  as  effectively  as  it  is  without 
requiring  more  advanced  methods. 


APPENDIX   A 

The   problem  is    stated   by   Den   Hartog    as    follows:      The    toy 
woodpecker   demonstrated   at    the    end   of   the    last    lecture    can   be 
idealized   as    shown    in    Figure   A-I.      The    sleeve    m   of    length    Si 
fits    loosely    around   the    fixed   stem.      The   birdie    (mass   M) 
tilts    the    sleeve    slightly  within    the    clearance,    so    that  point 
contact    occurs    at   A   and  B.      When   M  stands    still,    this    is    a 
self-locking    frictional   system.      When   M  moves    upwards,    the 
torque    of   the    sleeve    is    released   and   the    sleeve    slips    down    a 
bit.      Assume    the    sleeve   mass    m    (including    the    arm  to   the    left) 
to  have    its    center   of   gravity    in    the    center   of   the    vertical 
stem.       Let    the    friction    coefficient  be    f,    the    same    value    for 
motion    as    for   rest.       Let   x  be    the    downward   displacement    of 
the    sleeve    and   let   y  be    the    downward   displacement    of  M   rela- 
tive   to    the    sleeve   m,    measured   from  the   position    of  static 
equilibrium.      The    system   is    self-locking    at    rest    (y   =    0)    and 
X   then    is    constant   in    time.      For  y  negative    (upwards)    there 
may  be    slip. 

(a)  Discuss  the  static  equilibrium  problem  and  write  an 
expression  for  the  net  downward  force  on  m.  This  expression 
is  used  only  when  positive;  for  negative  values  the  friction 
reverses    sign.      Here    is    the   nonlinearity ! 

(b)  Set   up    the    two   differential   equations    of  motion    of 
the    system.      There   will  be    two   such   pairs:       one    good  while   m 
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(a) 


(b) 


Figure   A-I 
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is    slipping    (variables    x,  ,y,);    another   good  when   m   is    at 
rest    (variables    x^  ,79)  . 

(c)  As   you  have    seen    in    the    model    demonstration    at   the 
end   of   Lecture    23,    there    exists    a   steady,    periodic,    motion 
which    is    shown    in    Figure    A-II    in    its    general    characteristics. 
The   period   is    co   t   =    6    ,    divided   into    two  parts.      First   the 
sleeve    slips    from  A    (00   t   =    0)    until   B(u   t   =    <i>   )  ,   when    is 
seizes   by    friction.      Between    (})      <    (jo  t    <    6    ,    the    sleeve    m   is 
locked   on    the    stem   and   only    the   bird  M  moves    (the   y   shown 
above    is    positive,    i.e.,    the   bird   is    down,   binding   the    sleeve 
on    the    stem)  .       Since    in    general    6      is   not   equal    to    Ztt    (but 
somewhat    greater  probably)    the   natural    frequency    is   not    00^. 
Study    Figure   A-II    until   you   clearly   understand   its   physical 
implications . 

(d)  Now,   by   exact,    piecewise- linear  method,    write    the 
general    solution    of   the    four   equations    for   ^tjYi     (good   from 
A   to   B    above)    and   for   y^n'^?    ^S°°'^    from   B   to   C)  .      Since    all 
individual   equations    are    linear,    this    can  be    done   without 
great    difficulty;    the    solution    contains    many    (up    to   6    or    7) 
integration    constants    C, ,    C^ ,    ...    . 

(e)  Set   up    the   boundary    (continuity)    equations    at    the 
points    A,    B,    and   C    (w   t    =    0  ,    cf)    ,    e    ).      These   will    involve    the 
unknowns    C,  ,    C^  ,    ...,    ^^,    6^.      Check    and  polish   until   you 
have    as    many    condition   equations    as    there    are    unknowns. 
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So  =  t^nt 


Figure    A- I I 
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Solution 

(a)    IVhen    the    sleeve    is    just    ready   to   slip    downwards,    but 
does   not   move    yet,    it   must   be    in   static   equilibrium   under 
forces    as    shown    in    Figure   A-III. 
Moment    equilibrium   about    the    center   requires: 

(Mg   +   ky)a   =   N£ 
or 

N   =    (Mg   +   ky)a/ii. 


The   vertically    downward   force    on    the    sleeve    then    is 


(M+m)g   +   ky    -    2fN   =    (M+m) g  +   ky 


2fa 


(Mg+ky) 


Mg    1 


m 

M 


2fa 


I 


ky 


1    - 


2  fa] 


I 


(A-0) 


Ma  1 

The   toy's    dimensions    are    approximately  —  =    10;    j  =    4;    f   ~   j, 

and   it    is    clearly   self-locking  when  not    vibrating,    i.e.,    for 

y   =    0.      Then    the    downward   force   must   be   negative    or   at    least 

not  positive: 

^m_2fa_-,l      .7    =    . 9      <r, 

With  these  values ,  the  downward  driving  force  (a)  becomes 
-^  Mg  -  ky ,  which  means  that  it  drives  downwards  only  for 
negative,  i.e.,  upward  values  of  y.   This  is  physically  clear, 
because  for  upward  y,  the  normal  force  N  is  diminished  and 
the  friction  lock  is  broken. 

(b)  and  (c)  Now  consider  expression  (A-0)  to  be  the  down- 
ward driving  force  on  the  sleeve  during  slipping  motion 
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£^fM 


-^  N 


mg 


.ri£s^ 


il 


""S^ 


Figure   A-III 
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(x  >  0)  .  iNHienever  x  =  0,  the  expression  (a)  is  meaningless, 
because  the  friction  is  no  longer  f N ,  but  smaller  than  that. 
The    differential    equations    of  motion    are: 


for   X    >   0 
(from  A  to   B) 


for   X  =    0 
(from  B   to   C) 


mx^    =    Mg(l    +   J   -  2fa^    +   ky^(l    -    ^^)            (A-1) 

M(x-^   +   y^)    +   ky^  =    0                                                 (A-2) 

X2    =    0  (A- 3) 

My*      +   ky^    =    0  (A- 4) 


(d)    To    solve    these    differential   equations    (all    linear) , 

2    ^   k 
n  M^  n 


use    the   notation    oj^      =   rr.      This    co^    is    the   natural    frequency 


of   the   bird   on    its    spring  when   the    sleeve    is    fixed   and  not 
slipping;    it    is   not    the    frequency    of   the   periodic   motion 
including   slip.      Substitute    x,    from  equation    (A-1)    into 
equation    (A-2) ,    divided  by  M: 

••  .  2  ^  M  ,,  ^  m  2fa-^  ^  k  f^  2fa.  f. 
y,  +  (jo  yn+—  g(l  +  ir7-  —n — )  +  —  yT  (1  -  — „ — )  =  0 
-'I  n   -^  1        m  ^^  M  Si    ■'         m  -^1^  I    ' 

or 

Now    let 

^        m    "^^  Si    ^  ^    ' 

Ma  1 

a   constant    of   the    system.       (For   —  =10;-g-=4;f   =   -^we   have 

-A^    =    -9    or  A  =    +3)  .      Then: 
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with  the  solution: 

y^  =  C^e^'^^  -  C^e-^"^^    -  -\  .  (A-5) 

n 

To  find  X,,  we  return  to  equation  (A-2) : 

1       n    i         ^  " 

x^  =  -(a2  .  D.^^^C^e^^^  .  C^e-^'^nt^  ,  g 

x^  =  -  ^  [C^e^^n^  -  C^e-^'^^]  .  g|i  .  C3t  -  C^.      (A-6) 

Now  we   pass    to   equations    (A-3)    and    (A-4)    for   the    X2 ,72    func- 
tions   for   the   no-slip    domain.      We    leave    equation    (A-3) 
without    further  work    out,    and   integrate   equation    (A-4)    to: 


y-    =    C-    sin   w   t    +    C^    cos    m   t 
■^2  5  n  6  n 


(A-7) 


The  equations  (A-5),  (A-6),  (A-3),  and  (A-7)  are  the  answer 
to  part  (d) . 

(e)  Now  the  boundary  conditions.   Referring  to  the  figure 
on  page  46,  we  see  that  the  conditions  for  the  y  motion  are: 


^^l\   t=0    ^^I'^u}   t=e 
n  no 


(A-8) 


(yi).  t=o  =  ^^2^.  t=e  ^^-'^ 

n  no 

(^I    -)  =  fv  1  (A-10) 

no         no 
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no        no 
The  geometric  conditions  for  the  x  motion  are: 

n 

n    o 

There  is  one  more  (and  very  important)  condition  for  the  x 
motion  of  a  mechanical  nature.   Just  before  A  (before  co  t  =  0) 
the  sleeve  was  self- locking ,  but  at  point  A  it  just  ceases  to 
lock.   The  net  driving  force  at  A  then  just  breaks  away  from 
zero,  so  that 

n 
The  acceleration  just  before  B  (oo  t  =  *  )  is  not  necessarily 
zero.   Summarizing,  we  have  seven  condition  equations  (A-8) 
to  (A-14)  and  apparently  we  have  eight  unknowns: 

But  we  notice  that  the  constant  C.  only  occurs  in  equation 
(A-6)  for  the  x,  motion.   The  seven  condition  equations  (A-8) 
to  (A-14)  do  not  mention  x,  itself,  but  only  x,  and  x, ,  and 
we  see  from  equation  (A-6)  that  the  constant  C.  drops  out  in 
the  differentiation.   Thus,  we  have  seven  unknowns 

C^,C2  ,C^,C^  '"^G'^o'^o 

in  the  seven  equations  (A-8)  to  (A-14).   In  what  follows  the 

letters  ^      and  6   occur  many  times;  from  now  on  we  will 
o      o  ^ 


just   write    (p    and    6   without    subscripts,    for   simplicity.      Sub 
stitution    of   equations    (A-5) ,    (A-6) ,    (A-3),    and    (A-7)    into 
the    condition   equations    (A-8)    to    (A-14)    gives: 

C,    +    C„    -    -^  =    Cr    sin   0    +    C,    cos    6 
12^25  6 

n 
A(C^    -    C2)    =    C      cos    e    -    C,    sin   6 

C^e^*    +    C^e'^* ^  =    Cg    sin   <i>    +    C^    cos    <t> 

n 

A   C,e"^*    -    A   C„e'^*    =    C^    cos   (J)    -    C,    sin   4) 
i  Z  b  0 

-(a2   .   D^a  (Cje**    -    C,e-»)    .  S-*    .   C3  =   0 

n 

-(A^    +    l)\^(C^    +    C2)    +    g    =    0 

Eliminate  C^  between  (A-12)  and  (A-13) ,  leading  to: 

^-   [C^(l-e^*)  -  C2(l-e"^*)]  +  -^  c}>  =  0 .  (A-12*) 

CO 

n 


A^  +  1 


Solve    for   C^    and   C,    from    (A-8)    and  (A-0) : 

Cg    =    (C-j^    +    C2    -    -£2-)sin    e    +   A(C-j^    -  C2)cos    ( 

n 

S    "    '^'^1    "^    "^2    "    -^)cos    e    -   A(C^    -  C2)    sin 


(A- 

■8) 

(A- 

•9) 

(A- 

■10) 

(A- 

■11) 

(A- 

■12) 

(A- 

•13) 

(A- 

■14) 

n 
and   substitute   these    first    into    (A-10) ,    and    later   into    (A-11) 


C^e'^*+C2e'^*-    -^  =  ^'^l'^'^2"    "S"^  ^^^"    6    sin    ({)   +    cos    6    cos    (f)) 
w  00 

n  n 

+   A(C, -C2)  (cos  e    sin   <i>    -    sin   6    sin    cfi)  . 
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C^e^'^+C   e"^*-    -^  =    (C^  +  C2-    -^)  cos  (6- cj)) -A(C^ -C^)  sin  (0- ({))  . 

0)  "  CO 

"  ^  (A-10*) 

A  C^e^*-A  C   e"^*   =    ^^I'^'^z'    -^^  ^^^^    9    cos    cf)    -    cos    6    sine})) 


% 


+   A(C,-C2)(cos    e    COS    (f)   +    sin    6    sin    4)) 
A   C   e^*-A   C2e'^*   =    '^^l"'^2"    -S-)sin  (9- cj)) +A(C^ -C2)  cos  (e- 4>)  . 

CO 

"  (A-11*) 

By    this    time    the    set    of   seven   equations    (A-8)    to    (A-14)    have 
been  whittled   down   to    four   equations    (A-10*) ,    (A-11*) , 
(A-12*)    and    (A-14),    involving   the    four   variables    C,  ,    C2 ,    (J) 
and    (9-4)).      We   now   solve    for    C^    and   C^    from    (A-14)    and    (A-12*) 

r  g  l-Acj^-e-^'^ 

1  "    03   2(A^1)     '    2-e^*-e-^* 

n    ^ 

r      -  R  .    l-A»-e^^^ 

2  "    ^   2      2^,3         2-e^*-e-^* 

n 

with   the    combinations : 

C-    +    C-    =   — 9—^9 =   P       [see   equation    (A-14)] 

•■■  "^        w    ^(A^  +  1) 


2A*.e^*^e:^ 


Ci    =    4    =   P      2-e^*-e-^*         ' 


where   we   have    introduced   the    abbreviation 


n   =   g 

p  7        7 

CO    ^(A-^  +  l) 
n    ^  ^ 

Then  we  substitute  C,  and  C2  and  the  two  combinations 

(C  +C2)  and  (C^-C2)  first  into  equation  (A-10*)  and  later 

into  (A-11*): 
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=    [p    -    pCA'^  +  Dlcos    ^   -    Ap    "'    '^  L    "-A(i)      ^^"^    ^' 

2-e    ^-e 

where   we   have    used   the    abbreviation: 

i|j  =   e   -   4). 

Dividing  all  terms  by  p  and  multiplying  them  with  the  denomi- 
nator of  most  o£  the  terms  gives : 

e^*  -  Ac^e^*  -  1  -    e-^*  -  Ac^e"^*  =  1  =  (A^.l)  (2 -e^^-e"^*) 

=  -A2(2-e^*-e-^*)cos  ^  -  (-2A4,^Ae^*-Ae  ■^*)sin  ^ 

or 

2    cosh   A<i>    -    2A(t)   sinh   Acj)    -    2    -    2A^(l-cosh   H)    -    2(l-cosh   A^,) 

=    -2A^Cl-cosh   A(|))cos    \|j   +    2A(Act)-sinh   A(]))sin    ^ 
or ,    divided  by    2 : 
[-2    -    A^    +     CA^  +  2)cosh    A(j)    -    Ac})    sinh   AcJ)] 

=    [A^Ccosh   A(|)-l)]cos    ij;   +    [A(A(})-sinh   A({))]sin   i|j.  (A-10**) 

In   the    same   manner,   we    find  equation    (A-ll**) : 
[-A   sinh   Acj)    +    A^cf)    cosh   Acj)] 

=    -[A^Ccosh   Ac})-l)]sin   ij;   +    [A(Acj)-sinh   AcJ))]cos    i|j.         (A-H**) 

These  are  two  equations  in  the  unknowns  (j)  and  ^,    and  they 
have  the  form: 
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T   =    ot   cos    4^  +    B    sin    ^ 

5  =    -a  sin    i|J  +    3    cos    i) 

where  ct,  6,  y,  and  6  are  complicated  functions  of  c|).  We  now 
eliminate  4>  from  the  above  pair  by  squaring  each  one  of  them 
and  then    adding: 

Y^    +    6^    =    (a^    +    6^)(cos^iJj  +    sin^4i). 

Our  final  result  in  the  analytical  part  of  the  investigation 
thus    is  : 

a^    +    B^    -    Y^    -    6^    =    0, 

where 

2 
a  =   A    (cosh   A(f)-1) 

B    =   A(A(t)-sinh   A(j)) 

Y   =    (A^  +  2)(cosh   A(J)-1)    -    A())   sinh  Ac}) 

2 

6  =    -A  sinh   A<i>   +   A   (p    cosh   A<}). 

Further  analytical  work  does  not  seem  to  yield  any  simplifi- 
cation,   so    that   we   now   turn    to   numerical    computation. 

The   quantity   A,   which   must   be   positive,    is    a   character- 
istic parameter   of   the    system,    expressing   the    amount    of   dry 
friction    locking.       For   A   =    0    [see    equation    (A-0)]    the    fric- 
tion   f   is    so    small   that    the    sleeve    is    just    ready   to   slide 
down  with    zero   acceleration.      For  A  >    0,    there    is    a  self- 
locking   at    rest.      The    range    of  practical   values    of  A   appears 
to  be    from  0    to    3.       For   A  =    0,    all    four  expressions    a   to    6    go 
to      zero,    so    that   no    direct    conclusions    can   be    drawn    (we   have 
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to   employ   a    limit   process).      Hence,    take 
A  =    0.1,    1,    2 ,    3. 

Then   put    in    <i>   in   steps    of   one    degree    from   0    to    360    degrees 

2  2  2  7 

and   find   the    zeros    of    a     +    B      -    y      -    6    .       It    is   possible 

that    several    such    zeros   will   be    found.       Since    a   squaring 

operation   has    occurred    at    the    end   of   our   analysis,    not    all 

of   the    ({'-solutions   may   be    valid.       For   each    (|),    the   motion  has 

to   be    computed   and   inspected.      This    I   have   not    done,    and   if 

anyone   ever   does    it,    I    shall  be    grateful    if  he    or   she  will 

send  me    the    solution. 

This    concludes    Den   Hartog's    remarks.      There    are   probably 
several    reasons   why   the    analysis    is   not   more    fruitful.      With- 
out   considering    any    of   the    experimental    results    as    to    the 
relative    stability   of   the    two   models    in   Figure   A-I,   we    can 
see    in   the    light    of   Chapter    II    that   several   things    should  have 
been   included   in    the    above    analysis : 

1)  Losses    in   the    spring.       These    are    assuredly  present 
and   do    influence    the   motion.      When   a   larger   than  necessary 
displacement    is    used  to   start    the   motion,    most    of   the    excess 
energy   seems    to  be    dissipated   in   the    spring   as    the    sleeve    is 
not    observed   to   start    dropping    for   several    cycles. 

2)  The    impulsive   nature    of   the    motion  when   the    sleeve 
locks.       In   this    light   boundary    condition    (A-11)    is    incorrect 
and  probably   should   read 


^^1^0.    t=*      "    %    =    ^^2^0)    t=* 
n      ^o  n      ^o 
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where  v   is  related  to  the  drop  of  the  system  as  in  Chapter 
II. 

3)  Energy  balance  considerations.   These  would  relate 
directly  to  whether  stability  is  theoretically  possible  since 
if  energy  input  per  cycle  is  not  equal  to  energy  loss  per 
cycle  the  motion  will  either  delay  and  cease,  or  grow  uncon- 
trollably. 


APPENDIX  B 

Let  us  consider  a  situation  where  there  are  two  constant 
energy  impulses  delivered  each  period  and  let  the  losses  be 
due  to  linear  (velocity  dependent)  friction.   Let  the  loga- 
rithmic decrement  per  half  cycle  be  denoted  by  d,  and  let  the 
phase  velocity  be  y.   Just  after  an  impulse  let  the  phase 
velocity  be  y,.   A  half  cycle  later  it  will  be 

y   =  y-,e"^/2  _  (B-i) 

'  o   '1 

Just    after   this    an    impulse    occurs   which    changes    the   velocity: 

y/  -  y„^  =  h^  CB-2) 

7 
The  quantity  h   is  a  constant  representing  a  constant  energy 

input  during  the  impulses.   Just  after  the  impulse  the 

velocity  is  thus 


yo  =  [y,   +  h  ]  ' 


or 

y2  =  [Yi'e-^  -  h2]l/2  (B-3) 

using  (B-1).   We  require  that  y^  equal   y2  for  periodic 
motion  so  that  we  have 

yi-r2'y-  77-^172  •  ("-^^ 
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This  defines  a  closed  phase  path  corresponding  to  periodic 
oscillations.   The  motion  can  be  further  investigated  by  con- 
structing the  graph  of  the  sequence  function  (B-3)  [4,  p. 
433],  sometimes  called  "Lamerery's  diagram"  [2,  p.  161].   This 
is  shown  in  Figure  B-I.   The  vertical  axis  plots  values  of 
y^  and  the  horizontal  axis,  y^.   For  a  stationary  oscilla- 
tion y-.  equals  y^  and  this  is  represented  by  a  straight  line 
through  the  origin.   If  the  oscillations  are  not  stationary 
the  sequence  of  amplitudes  can  be  determined  by  (B-3)  and 
Figure  B-I.   If  the  initial  amplitude  is  yi ,  then  the  subse- 
quent amplitude  will  be  yl ,  and  the  corrected  amplitude  will 
be  Yj    on  the  y-,  axis.   The  next  amplitude  will  be  yl   on  the 
y^  axis  and  the  sequence  of  amplitudes  will  continue  until  y 
is  reached,  at  which  point  a  stable  limit  cycle  exists. 
Physically  this  corresponds  to  a  situation  where  initially 
the  energy  introduced  is  greater  than  that  dissipated  and  the 
oscillations  grow  in  amplitude  until  a  stable  energy  balance 
is  reached. 

The  important  point  to  note  is  that  y^  can  be  made 
arbitrarily  small  and  the  oscillations  will  still  grow  until 
a  stable  limit  cycle  is  reached.   This  implies  soft  excitation 
is  possible,  a  phenomenon  which  is  not  observed  physically. 
In  fact,  Andronov  concludes  his  analysis  [2,  p.  172]  with  the 
statement,  "Again,  a  model  with  linear  friction  does  not 
explain  the  need  for  an  initial  finite  impulse  to  start  the 
clock  and  must  therefore  be  rejected.   It  is  necessary  to 
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>? 


n^Cy.^e-^^v.^]^ 


X 


Figure    B-I 
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assume  that  the  clock  is  a  self -oscillating  system  with  dry 
friction. " 

The  fact  that  Coulomb  friction  is  appropriate  can  be 
seen  by  considering  a  phase  trajectory  of  the  actual  wood- 
pecker.  This  is  shown  in  Figure  B-II.   This  trajectory 
terminates  on  the  x  axis  between  the  points  (-b  -f  ,0)  and 

fb  +f  ,0)  which  characterizes  the  zone  of  static  friction. 
*-  o   o '  -^ 

When  a  trajectory  does  not  emerge  from  this  zone  the  motion 
ceases.   Physically,  this  corresponds  to  the  situation  where 
the  amplitude  of  the  bird  mass  is  not  sufficient  to  cause  the 
sleeve  to  unlock  initially.   The  quantity  b   used  in  Chapter 
II  can  now  be  seen  as  a  parameter  describing  the  friction 
lock  of  the  sleeve  as  far  as  its  effect  on  the  motion  is  con- 
cerned. 


Y 


Figure    B-II 


APPENDIX  C 

The  following  are  examples  of  data  collected  from  the 
videotape  made  of  the  woodpecker  motion.   The  times  listed 
were  read  from  an  electronic  timer  with  a  digital  readout 
and  are  accurate  to  ±0.005  second.   The  hundredths  position 
was  sometimes  difficult  to  read.   The  elapsed  time  is  more 
definite  with  each  reading  occurring  one-sixtieth  of  a  second 
after  the  last.   The  position  readings  refer  to  positions  in 
centimeters  ±0.15  read  from  a  labeled  grid  behind  the  wood- 
pecker.  Vertical  displacements  are  labeled  y  and  horizontal 
ones,  X.   The  first  series  of  readings  are  for  a  bird  mass  of 
11. 30  grams . 

Sleeve  y    Bird  x    Bird  y    Time 


0.2 

7.2 

8.8 

3.80 

0.2 

7.9 

8.2 

3.82 

0.2 

8.9 

7.4 

3.84 

0.2 

10.0 

5.9 

3.86 

0.2 

10.6 

4.8 

3.87 

0.2 

11.1 

3.2 

3.89 

0.3 

11.3 

3.0 

3.91 

0.4 

11.3 

1.0 

3.92 

0.5 

11.2 

0.1 

3.94 

0.5 

10.0 

-0.7 

3.96 

0.5 

11.0 

-0.7 

3.97 

0.5 

11.0 

-0.2 

3.99 

0.5 

11.2 

0.2 

4.01 

0.5 

11.3 

1.6 

4.02 

0.5 

11.1 

3.1 

4.04 

0.6 

10.7 

4.9 

4.06 

0.6 

9.8 

6.3 

4.07 

(continued) 
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Sleeve  y 

Bird  X 

Bird  y 

Time 

0.6 

8.  7 

7.8 

4.09 

0.6 

8.0 

8.5 

4.11 

0.6 

6.9 

9.1 

4.12 

0.6 

6.8 

9.2 

4.14 

0.6 

6.1 

9.1 

4.16 

0.7 

8.0 

8.5 

4.17 

0.7 

9.0 

7.6 

4.19 

0.7 

9.9 

6.1 

4.21 

0.7 

10.8 

4.9 

4.22 

0.7 

11.0 

3.5 

4.24 

0.7 

11.2 

2.3 

4.26 

0.8 

11.2 

1.2 

4.27 

0.85 

11.0 

0.3 

4.29 

0.85 

10.9 

-0.1 

4.31 

0.  85 

10.8 

-0.2 

4.32 

0.85 

11.0 

0.1 

4.34 

0.85 

11.2 

0.9 

4.36 

0.85 

11.4 

2.0 

4.37 

0.90 

10.1 

3.5 

4.39 

1.00 

10.9 

5.2 

4.41 

1.00 

10.1 

7.0 

4.42 

1.00 

8.9 

8.2 

4.44 

1.00 

7.9 

9.0 

4.46 

1.00 

7.1 

9.8 

4.47 

1.00 

6.8 

9.8 

4.49 

1.00 

7.1 

9.7 

4.51 

1.00 

7.9 

9.0 

4.52 

1.00 

8.8 

8.2 

4.54 

1.00 

10.0 

6.9 

4.56 

1.00 

10.8 

5.2 

4.57 

1.00 

11.1 

4.0 

4.59 

1.00 

11.6 

2.9 

4.61 

1.10 

11.5 

1.8 

4.63 

The    following   data  were    taken    for   a  bird  mass    of   11.30 
grams   with   the   sleeve    locked   in  place. 

Bird   X  Bird  y  Time 

12.1  7.3  2.08 

12.0  6.5  2.10 

12.0  6.0  2.11 

11.9  5.9  2.13 

12.0  6.2  2.15 

12.0  6.8  2.17 

12.0  7.7  2.18 

(continued) 
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Bird  X    Bird  y    Time 


11.7 

9.0 

2.20 

11.1 

10.1 

2.21 

10.5 

11.1 

2.23 

9.8 

12.0 

2.25 

9.0 

12.9 

2.27 

8.  7 

13.2 

2.28 

8.6 

13.2 

2.30 

8.9 

13.0 

2.31 

9.6 

12.4 

2.33 

10.2 

11.4 

2.35 

11.0 

10.4 

2.37 

11.4 

9.4 

2.38 

11.9 

8.3 

2.40 

12.0 

7.3 

2.42 

12.0 

6.9 

2.43 

12.0 

6.4 

2.45 

11.9 

6.1 

2.47 

12.0 

6.3 

2.48 

12.0 

6.9 

2.50 

Figure  C-I  shows  a  graph  of  some  of  the  data  points  for 
the  position  of  the  sleeve  versus  time  and  Figure  C-II  shows 
a  graph  of  the  speed  of  the  sleeve  versus  time  obtained  by 
calculating  the  slope  of  the  line  in  Figure  C-I. 
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